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Discovery of a quiescent neutron star binary in the globular 

cluster M13 
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Abstract. We have discovered with XMM-Newton an X-ray source in the core of the globular cluster M13, whose 
X-ray spectral properties suggest that it is a quiescent neutron star X-ray binary. The spectrum can be well fitted 
with a pure hydrogen atmosphere model, with Too = 76 ± 3 eV, Roo = 12.8 ± 0.4 km and an X-ray luminosity of 
7.3 ± 0.6 x 10 32 ergs s" 1 . In the light of this result, we have discovered a strong correlation between the stellar 
encounter rate and the number of quiescent neutron stars found in the ten globular clusters observed so far by 
either XMM-Newton or Chandra. This result lends strong support to the idea that these systems are primarily 
produced by stellar encounters in the core of globular clusters. 
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1. Introduction 

Globular clusters (GCs) are known to harbour both bright 
and dim X-ray sources. Bright X-ray sources with lumi- 
nosities above ~ 10 36 ergs s _1 are commonly agreed to be 
neutron star Low-Mass X-ray Binaries. There are twelve 
such sources known in GCs; 11 of them have shown type 
I X-ray bursts, the unmistakable signature of an accreting 
neutron star. On the other hand, dim GC X-ray sources 
are much more numerous (for example, dozens of such ob- 
jects exist in to Cen, see iGendre et al~ll2003|) . but their 
nature is still debated. They have maximum luminosi- 
ties of ~ 10 33 ergs s _1 and there is growing evidence 
that they are a variety of different objects. Some of the 
fainter ones have been identified as Cataclysmic Variables 
or active binaries (RS CVn, BY Dra). Others, have been 
associated with radi o millisecond pulsars (MSPs) (e.g. 
iGrindlav et al. Ill99l|) . Some of the 'brighter' of the dim 
X-ray sources have recently been proposed to be quiescent 
neutron star X-ray binaries (qNSs) on the basis of their 
X-ray spectral properties (47 Tuc, Grindlay et al. 2001; u> 
Cen, Rutledge et al. 2002, Gendre et al. 2002; NGC 6440, 
Pooley et al. 2002b). 

Outside GCs, qNSs are found to have sof t X-ray spec- 
tra a nd luminosities up to 10 33 ergs s _1 (e.g. IKong et al. I 



2002). During ou tbursts, reactions, d eep in the crust, heat 
the neutron star ijBrown et al. Ill998^ . Between outbursts, 
the heated surface radiates a thermal spectrum, emitted 
by a neutron star hydrogen atmosphere (NSA) . The ther- 
mal component of the quiescent X-ray spectra of two neu- 
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tron star X-ray binaries which are not in a globular cluster 
(Aql X-l and Cen X-4) are well fitted by NSA models, pro- 
viding strong support for this theory. In contrast, quies- 
cent black hole transients have hard power law like spectra 
and low luminosities (down to 10 30 ergs s" 1 ). Accretion 
via an advection dominated flow is thought to be responsi- 
ble for the observed hard an d weak emission ijKong et al. I 
l2002tlHame"ur7 et al. ll2003l and references therein) . Thus 
the combination of a soft thermal X-ray spectrum and a 
luminosity above ~ 10 32 ergs s _1 has been used to claim 
the detection of qNSs in GCs. Th e existence of such sys - 
tems in GCs was first proposed bv lVerbunt et al."l l)l984|) . 
It is generally assumed that they are formed through stel- 
lar encounters (tidal capt ure or exchange encounters) in 
the dense cores of GCs fsee lHut et al. Il992t for a review). 

ROSAT HR I and PSPC observations l|Fox et al. Il996t 
IVerbunt I Eiol . revealed several dim X-ray sources in 
M13. The HRI observations showed two core sources, 
whose unresolved PSPC spectrum could be fitted with 
a 0.9 keV thermal bremsstrahlung model (L0.5_2.5keV ~ 
2.5 x 10 32 ergs s _1 ). M 13 is somewhat remarkable as it con- 
tains also four MSPs l|Tavlor et al. lll993T: iRansom et al. I 
120021) . 



2. Observations and results 

The globular cluster M13 was observed by the XMM- 
Newton EPIC cameras on 2002 January 28 and 30, us- 
ing the Full Frame Window and a medium filter. The 
total length of the observation was ~ 37 ks. We ana- 
lyzed the data using the XMM-Newton Science Analysis 



2 



Gendre, Barret, Webb: XMM-Newton observation of M13 









- 


I / 




N 




/ 

/ 




^ \ 


\ 

---\ 

1 

/ ; 




\ ^ 


['•J 

f ' — 

\ ° ♦ 






\ 

\ 


\ * i / 




\ 


■ ; / 


/ - 


c 


\ 








, , , ; , p. , ; , , , 


- 



16" 41"" 50° 



RIGHT ASCENSION 



Fig. 1. A contour image of the center of the field of view 
of M 13. The core and half mass radii are shown as solid 
and dashed lines, respectively. Pr evious id entifications are: 
ROSAT sources (filled circles. IVerbunt 11200 ilk faint UV 
sources (filled s quares, iFerraro et al. Il997lk radio objects 
(open squares, I.Tohnston et al. I Il99llk and MSPs with 
known positions (filled diamonds . iTavlor et al. Il993i) (the 
most central MSP has positional uncertainties larger than 
the image displayed. All other objects have uncertainties 
smaller than the symbol size). 



Software (SAS) version 5.3.3. Initially we considered the 
two segments of the observation separately. We used the 
SAS tasks emchain and epchain to calibrate the raw data, 
flag bad pixels and filter for non-astrophysical events. The 
background was found to be variable and relatively high. 
Removing the periods of unstable background left 9 and 8 
ks respectively. The filtered event files were merged using 
merge and then we extracted images and spectra. 

The source detec tion procedure used is described in 
iGendre et al. I l|2003l ) . Briefly, it combines a wavelet detec- 
tion algorithm applied to a 0.5 - 5.0 keV band image, and 
a maximum likelihood fitting of the source candidates. A 
conservative maximum likelihood threshold of 12 was cho- 
sen. The absorption-corrected limiting flux is 3.6 x 10~ 15 
ergs s _1 cm' 2 corresponding to a luminosity of 2.6 x 10 31 
ergs s -1 for a source located at the center of the field of 
view having a 0.6 keV blackbody spectrum. 

We found an extended source within the core radius. 
The extension appeared to be due to two sources that were 
not quite resolved. The MOS cameras have the advantage 
of having a smaller pixel size than the PN camera: 1.1" 
for MOS llTurner et al. Il200l|) as oppose to 4.1" for PN 
ijStriider et al. Il200l|k When processing the event file, the 
task emevents converts the event position (RAWX and 
RAWY) into camera coordinates in units of 0.05". This 



step includes randomization (within a CCD pixel) of the 
event to avoid the Moire effect. By default, the xmmse- 
lect task, assumes a binning factor of 87 to produce MOS 
images with the same pixel size as PN images. We have 
produced MOS images with a binning size of 20, corre- 
sponding to a pixel size of 1". The image smoothed with 
a simple 2D-Gaussian of a = 2.0" is shown in Fig^ 

To summarize, we detected 2, 5 and 77 sources within 
the core radius, half mass radius and field of view re- 
spectively. Using t he Log N-Log S c urve of extragalactic 
sources reported in lHasinger et al. 1 l|200l|) and the limit- 
ing flux in the three regions given above, we determined 
the expected number of background sources to be 0, 1 and 
72, respectively. The positions and errors of sources found 
within twice the half-mass radius are given in Tabled It 
goes beyond the scope of this paper to investigate the er- 
ror box content of each of the XMM-Newton sources. It is 
however worth mentioning that only 10 of the 11 ROSAT 
sources that should have been detected by XMM-Newton, 
given the limiting flux and field of view of the observation, 
were detected. We failed to detect the ROSAT source Gb. 
This source must have therefore varied in flux by at least 
a factor 10 between the ROSAT and XMM-Newton ob- 
servations. 

Table 1. XMM-Newton sources detected within twice 
the half mass radius (HMR) (CR = Core Radius). The 
positional error includes the statistical error (90% confi- 
denc e, estimated from sr ematch) and a systematic error 
of 4" ll.Ta,r l sen et a ,l.l l200llk 



R.A. 

h m s 



Dec. 

O / // 



Error Location 



16 41 37.9 


36 28 26.3 


7.94 


HMR 


16 41 42.7 


36 28 06.7 


4.60 


CR 


16 41 43.8 


36 27 58.6 


4.27 


CR 


16 41 46.8 


36 27 29.5 


5.93 


HMR 


16 41 49.0 


36 26 44.6 


7.53 


2xHMR 


16 41 38.3 


36 26 27.0 


4.50 


HMR 


16 41 49.9 


36 26 18.7 


7.38 


2xHMR 


16 41 42.5 


36 25 53.0 


6.83 


2xHMR 



The two core sources are separ ated by only 15 " (see 
FigHJ. One is the ROSAT source Ga l|Verbunt l200llk The 
other should have been detected by ROSAT, if its flux had 
not varied. We reanalyzed the HRI data and determined 
that the source flux must have varied by a factor ~ 2 
between the ROSAT and the XMM-Newton observations. 

Two sources lying so close together complicates the 
spectral analysis. Normally, spectra are accumulated over 
a region of radius of 0.7' to include 85% of the source pho- 
tons. Such an extraction region is 3 times larger than the 
source separation. Using a radius of 0.7', we have extracted 
the spectrum of the brighter source by masking out a re- 
gion of radius 15", offset by 5" from the fainter source. We 
have estimated that only 6% of the counts in the spectrum 
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of the brighter source come from the fainter source, insuf- 
ficient to affect the results of our spectral analysis. Spectra 
were extracted from the EPIC-PN and the two MOS cam- 
eras. We binned these spectra to contain at least 20 net 
counts per bin and generated ARF and RMF files with the 
SAS tasks arjgen and rmfgen. The spectra are extremely 
soft, with ~ 90% of the counts below 2 keV. We have tried 
to fit the combined spectrum with different single com- 
ponent models (blackbody and thermal bremsstrahlung) . 
Absorption by the interstellar medium was included in the 
fit but was found to be consistent with the expected value 
from th e optical extinctio n towards the cluster (1.1 x 10 20 
cm~ 2 , iDioreovski Ill993|) . We also tried t o fit the spec- 
trum with a pure hydro gen NSA model l|Pavlov et al. I 
Il992t IZavlin et ai~lll996[) . This model provides the best 
fit to our data. Assuming a neutron star mass of 1.4Mq, 
we derived T^ = 76t 3 3 eV and = 12.8i° * km, with a 
xl=0.55 (15 degrees of freedom (dof)). These parameters 
are similar to those determined for the proposed qNS in 
lo Cen (Rutledge et al. 2002, Gendre et al. 2003). 

We have retrieved and reanalyzed the ROSAT PSPC 
archival observations of M13 to determine whether the 
PSPC spectrum of source Ga could be fitted with the same 
model. A fit of the combined PSPC and XMM-Newton 
data revealed parameters consistent with those derived 
from the XMM-Newton data alone (xl = 0.90, 26 dof). 
We present the unfolded combined EPIC-PN and ROSAT 
spectrum in Fig. |21 The luminosity derived for the NSA 
model is (7.3± 0.6) x 10 32 erg s" 1 (0.1-5.0 keV), using the 
distance of 7.7 kpc (Harris 1999). 



The quiescent neutron star candidate 




0.1 1.0 
Energy (keV) 



Fig. 2. The unfolded EPIC-PN (blue points) and ROSAT- 
PSPC (red points) spectra of the qNS candidate and the 
best fit NSA model, one of the best qNS spectra to date. 



3. Discussion 

The most likely interpretation of the nature of the softest 
source in M13 is that it is a quiescent neutron star low- 
mass X-ray binary. This interpretation is supported by the 
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Fig. 3. The number of qNSs versus collision rate (po' 5 f 2 ), 
shown with a linear fit n q NS ~ 0.04 x p\- 5 r 2 + 0.2. The 
collision rates have been normalized to 100 for NGC 6440. 



luminosity of the source, the softness of its X-ray spectrum 
and the fact that a NSA model yields a good fit and plau- 
sible parameters (radius and temperature) for the neutron 
star. A reasonable question to ask is whether we expect 
such a system in M13. In the disk, qNSs have a minimum 
X-ray luminosity of ~ 10 32 ergs s -1 l|Naravan et aT. 2002). 
As it has been already emphasized, if the same luminos- 
ity threshold also applies to globular cluster qNSs, the 
luminosity limit of XMM-Newton and Chandra observa- 
tions (typically around 10 30 — 10 31 ergs s _1 at the cluster 
distances) allows one to detect all the qNSs present in 
globular clusters. In table |3 we list the globular clusters 
already observed by either Chandra or XMM-Newton, to- 
gether with the number of qNSs reported in the literature. 

In globular clusters, the number of qNSs is expected to 
scale with the collision rate which is proportional to Po' 5 ^ 
for virialized clusters, where po is the central de nsity of 
the cluster and r c its core rad ius llVerbunt 112002^ . These 
values taken from the lHarris I l)l999|) catalog are listed in 
Table |21 In Fig [3J we plot the number of qNSs as a func- 
tion of the collision rate, normalized so that the value for 
NGC 6440 is 100. There is a striking correlation between 
the number of qNSs and the collision rate. The presence 
of one qNS in M13 is therefore not really surprising given 
that its collision rate puts the cluster in a region where 
one might expect either one or zero qNS. This remarkable 
correlation extends over more than 2 orders of magnitude 
and includes both core-collapsed and non core-collapsed 
clusters. This strongly supports the idea that qNSs are 
indeed prim arily produced by stellar encounters in globu- 
lar clusters l|Verbunt Il200l| . 

With the results of the observations reported here and 
the four already known MSPs, there are at least 5 neutron 
star systems in M13. This makes M13 the cluster with 
the fourth highest number of known neutron star systems. 
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Table 2. qNSs in GCs detected by Chandra or XMM- 
Newton. Parameters for the clusters are taken from the 
lHarris I l)l999(l catalog; the distance is given in kpc, the 
core radius in '. We indicate the log of the central density 
in units of L Q pc~ 3 . 



Cluster 


Distance 


Core 


Central 


qNSs 




(kpc) 


radius 


density 




47 Tuc 


4.5 


0.44 


4.77 


2 


u> Cen 


5.3 


2.58 


3.12 


1 


M 13 


7.7 


0.78 


3.33 


1 


NGC 6366 


3.6 


1.83 


2.42 





NGC 6397 


2.3 


0.05 


5.68 


1 


NGC 6440 


8.4 


0.13 


5.28 


4-5 


M 28 


5.7 


0.24 


4.75 


1 


M 22 


3.2 


1.42 


3.64 





NGC 6752 


4.0 


0.17 


4.91 






References, top to bo tt om: iGrindlav et al. I |2001aL l2002jh 
iRutledee et al. I J2002I) . iGendre et al. I (120031); this work; 
Webb, B arret, Gendre , in pre pa ration; |&ijid la^ T et_alJ 



J2001bl): iPoolev et al. 1 (l2002bl): 


Becker et al. 1 J2003D: 


IWebb et al. 1 (120021): andlPoolev et al. 


J20024. 



The retention of such a large number of neutron stars in a 
cluster with a relatively low central density remains to be 
explained (see for a compreh ensive study of neu tron star 
retention in globular clusters. IPfahl et al. lEoO^ . 

4. Conclusion 

We have reported the likely discovery of a quiescent neu- 
tron star in M13. We have also shown for the first time 
that there is a strong correlation between the stellar col- 
lision rate and the number of qNSs in the 9 globular clus- 
ters observed by either XMM-Newton or Chandra. This 
lends strong support to the idea that these systems are 
primarily produced by stellar encounters in the core of 
globular clusters. More observations, in particular with 
XMM-Newton are being planned and should enable us to 
test the strength of this correlation. 
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